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Abstract 
Calculations  with  a  synoptic  case  study  show  that  long-wave 
radiative  cooling  tends  to  reduce  the  available  potential  energy, 
especially  in  the  upper  troposphere.    Synoptic-scale  precipitation 
amounts  resulting  from  destabilization  of  clouds  by  long-wave  cooling 
are  computed.      These  range  up  to  1.4  mm  in  12  hours.    This  de- 
stabilizing effect  may  be  important  in  explaining  the  nocturnal  maximum 
of  precipitation  over  the  sea.    It  may  also  contribute  significantly  to 
cyclone  development. 


1 .      Introduction 

This  paper  is  a  continuation  of  Scientific  Report  No.   1  published 
earlier  under  this  Research  Project  (Danard  1968). 

Calculations  of  Katayama  (1967)  show  that  long-wave  radiation, 
absorption  of  solar  radiation  and  release  of  latent  heat  are  about  equally 
important  in  the  annual  mean  generation  of  available  potential  energy 
over  the  northern  hemisphere.     Section  2  is  concerned  with  the  first 
of  these  processes.    Cloud-top  cooling  reduces  temperature  contrasts 
if  the  cloud  mass  is  in  the  warm  air. 

The  existence  of  a  nocturnal  maximum  in  frequency  of  precipitation 
over  the  sea  has  been  demonstrated  by  Kraus  (1963).    He  studied  the 
records  of  weather  ships  in  the  eastern  Pacific  and  North  Atlantic 
Oceans.    With   8-10  years  of  data  for  July,  he  found  that  precipitation 
was  reported  in  19  per  cent  of  all  nocturnal  synoptic  reports  compared 
to  12  per  cent  of  the  daytime  reports. 

Kraus  proposed  that  during  the  day,  absorption  of  short-wave 
radiation  reduces  the  production  of  liquid  water  in  clouds,  thereby 
lowering  the  probability  of  precipitation.    At  night,  this  effect  is  ab- 
sent. 

The  nocturnal  maximum  in  precipitation  occurrence  over  the  central 
United  States  in  summer  was  studied  by  Bleeker  and  Andre  (1951). 
Their  explanation  was  that  low-level  differential  cooling  at  night  due 
to  variations  in  topography  would  result  in  upward  motion  between  the 


Appalachians  and  the  Rockies. 

In  Section  3  it  is  suggested  that  an  important  factor  is  the  de- 
stabilizing influence  of  long-wave  cooling  at  night  from  cloud-top.     The 
amount  of  precipitation  associated  with  convective  overturning  is  com- 
puted.   The  destabilization  may  also  augment  the  large-scale  vertical 
velocity.    This  would  further  increase  the  precipitation  as  well  as 
accelerate  cyclone  development. 


2  .    Effects  of  long-wave  radiation  on  the  available  potential  energy. 
Consider  the  thermodynamic  energy  equation 

Sa  H     r,  RH  ,,v 

F-V.7«+ffM  +  —  (1) 

P 

where    cr  =  -a9  j>n  Q/dp,     a.1  =  dp/dt,    and    H    is  the  rate  of  input  of 
heat  per  unit  mass .    In  (1),  let   a  =  a  +  a',  o=a  +  a',  etc. ,  where 
the  bar  denotes  the  average  over  an  isobaric  surface  of  area    S    and 
the  prime  the  departure  therefrom.    Eq .   (1)  may  then  be  transformed  to 
(Danard  1966  or  Wiin- Nielsen  1964) 


V-|A       J-  jPbj  v.7a.2dS..dp+I   S 


J    V  •  7  a'^dS'-dp  +-    J       J    a'to'dSdp 


where 


1     nPb  r    a'CT1  id'  ,10  .,  R        rPb  r    a'  H'    ,c   . 

+  -  — = —    dS  dp  +  — —  — dS  dp 

q    c  J  Q  Cg°  Jcrp 

Pt  s  p y      P.  s 

(2) 


1  Ph  rv'2 

A"£    I  *=T    dSdP  (3) 

g       p  s      a 
yt 


will  be  referred  to  as  the  available  potential  energy  of  the  region  (c  f , 
Eq  (10)  of  Lorenz  (1955)  or  Eq  (2.2)  of  Wiin-Nielsen  (1964)).    In  (2) 
and  (3),  p    and    p      denote  the  pressures  at  the  top  and  bottom  of  the 
region  of  interest. 

If  a  quasi-geostrophic  balance  exists,  co1  will  be  related  to    H' . 
However,  attention  will  be  restricted  here  simply  to  the  final  term 


in  (2),  representing  the  generation  and  designated  as  G: 

p 

G  =  —     J        J      -3—  dS  dp  (4) 

P  Pt      s       cr  p 

In  (4),  H'  will  be  identified  with  long-wave  radiative  effects  com- 
puted in  Section  6  of  Sci .  Rep.  No.  1.    The  area    S    is  bounded  by  the 
latitudes  28  .  7N  and  51 .  3N  and  by  the  longitudes  70W  and  115W.     Com- 
putations are  performed  for  three  layers:    900-700  mb,  700-500  mb,  and 
500-300  mb .    Results  are  given  in  Table  1 . 

In  the  lower  layers,    G    has  small,  positive  values.    This  is  due  in 
large  part  to  the  absence  of  cooling  in  the  lower  portion  of  the  main 
cloud  mass  (see  Figs.  8  and  9  of  Sci.    Rep.  No.  1.)    In  the  lowest 
layer,  there  is  also  the  effect  of  cooling  from  the  top  of  stratocumulus 
cloud  in  the  cold  air  over  Kansas  and  Oklahoma  (see  Fig.  9  of  Sci.  Rep. 
No.  1).    However,  in  the  uppermost  layer,    G    is  negative  and  fairly 
large  in  magnitude.    This  is  due  to  radiation  from  cloud-top  in  the  warm 
air.    Thus  long -wave  radiation  tends  to  reduce  the  temperature  contrast 
in  the  upper  troposphere,  a  not  unexpected  result.    Assuming    G    con- 
stant and  neglecting  the  other  terms  on  the  right  side  of  (2),  the  quo- 
tient A/ JGj  represents  the  time  required  to  effect  a  change  in   A 
equal  in  magnitude  to  the  initial  value.    Whereas  this  time  is  of  the 
order  of  weeks  or  longer  in  the  lower  layers,  it  is  of  the  order  of  days 
in  the  highest  layer.      This  suggests  that  long  wave  radiation  from  cloud- 
top  may  be  important  in  the  upper  troposphere  even  for  short  time  periods . 
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Table  1 .    Comparison  of  terms    A    and    G     .     See  Eqs  (2)  -  (4) . 


Term  Units  900-700  mb        700-500        500-300      900-300 


A 

in25 

10       ergs 

6.04 

7.36 

4.97 

18.4 

G 

m20                    -1 
10       ergs  sec 

0.26 

0.08 

-1.66 

-1.32 

a/|g| 

days 

27 

104 

3.5 
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3 .      Effects  of  long-wave  radiation  from  cloud-top  on  precipitation  rates 

The  basic  principles  in  computing  the  convective  precipitation  are 
similar  to  those  employed  by  Manabe  et  al  (1965).    Consider  an  unstable 
atmosphere  initially  at  rest.    If  the  motions  arising  from  convective 
overturning  are  dissipated  as  heat,  it  may  be  shown  that 

6  J  C  T  dm  =  J6Q  dm  (5) 

where    6    refers  to  a  change  over  a  finite  time  interval,    dm    is  an 
element  of  mass  and    6Q    is  the  heat  added  per  unit  mass  from  all 
sources  other  than  friction.    In  (5),  the  integration  extends  over  the 
entire  atmosphere.     Eq  (5)  is  derived  in  the  Appendix.    In  a  dry  adiabatic 
transformation  the  sum  of  the  potential  and  internal  energies  (represented 
by  the  integral  on  the  left  side  of  (5))  is  conserved. 

Eq  (5)  is  applied  to  a  layer  cloud  extending  from  1000  to  400  mb. 
The  lapse  rate  is  assumed  initially  pseudo-adiabatic.    To  simulate  con- 
ditions during  the  night,  the  cloud  top  radiates  for  12  hours.    As  this 
proceeds,  the  cloud  becomes  unstable.    The  convective  overturning 
then  produces  precipitation  and  returns  the  lapse  rate  to  the  pseudo- 
adiabatic  value.    The  liquid  water  contents  of  the  initial  and  final 
states  are  assumed  equal  (i.e. ,  all  water  which  condenses  during  this 
process  falls  out  as  precipitation). 

Applied  to  the  above  problem  (5)  may  be  written 

Cd        P°  L      Po 
L   j        6Tdp    =   AQ+-    J'      6rdp  (6) 


g      °  p4  9 


P4 


where    p    =  400  mb,    p    =  1000  mb,    6T(<0)    and    6r(<0)arethe 
differences  between  final  and  initial  temperatures  and  mixing  ratios, 
AQ  (>0)  is  the  heat  lost  by  radiation,  and    L    is  the  latent  heat  of 
vaporization.    The  total  precipitation  is 

P  =  -JP°6r^  (7) 

J  g 

P4 

The  interval  between    p,    and    p      is  divided  into  six  layers  each 

4  o 

100  mb  thick.    The  integrals  are  then  expressed  as  sums  of  contributions 
from  these  layers.    With  a  known  initial  state,  Eq  (6)  is  thus  an  equa- 
tion in  12  unknowns,  the  final  values  of    T    and    r    at  the  mid-points 
of  the  six  layers.    Eleven  more  equations  are  needed.    Five  are  obtained 
from  the  final  condition  that 

se 

W  -  °  <8) 

where 

P 

is  the  equivalent  potential  temperature.    The  remaining  six  equations  are 

obtained  from  the  Clausius-Clapeyron  equation.    If    r. ,  e.    and    T. 

ill 

denote  initial  mixing  ratio,  vapor  pressure  and  temperature,  one  obtains 

0.622  r.L 

6r  =  ~^-     —2—    6T  (10) 

RT.2  P'ei 

i 

Eq  (10)  is  applied  at  the  six  levels. 


In  computing  the  radiation  lost  from  cloud-top,  equations  similar 

to  (2)  and  (13)  of  Sci.  Rep.  No.  1    were  derived.    The  level    p      was 

taken  as  400  mb  (cloud-top)  and     q      was  set  equal  to  the  saturation 

u 

value  corresponding  to  the  temperature  there.    The  upward  flux  at  cloud- 
top  is  thus 

F     =F.      (1  -  €   )  (11) 

u        bu  u 

where  F,        is  the  black  body  flux  corresponding  to  cloud-top  temperature, 

and    €     is  the  water  vapor  emissivity  for  the  layer  above  the  cloud.    The 
u 

heat  lost  by  radiation  is  ,  per  unit  area  , 

AQ  =  F    At  (12) 

where  A  t  =  12  hrs . 

The  method  described  above  is  applied  to  three  cases  where  the 

initial  state  is  saturated  with  a  constant  wet-bulb  potential  temperature 

8       of  0,   10  and  20  C.    These  would  represent  conditions  over  polar, 
w 

mid-latitude  and  sub-tropical  oceans,  respectively.    Results  are  shown 

in  Table  2. 

The  values  of    AQ    decrease     with  increasing  initial    8    .    This 

w 

is  because  the  decrease  in    ( 1  -  €   )  in  Eq  (11)  over-compensates  the 

increase  in  F,     .    Release  of  latent  heat  tends  to  offset  the  cooling  due 

to  long-wave  radiation.    This  results  in  fairly  small  values  of    6  8    . 

w 

Precipitation  amounts  in  Table  2  are  to  be  regarded  as  averages  over 

large  areas,  and  not  as  values  associated  with  individual  convective 
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Table  2  .    Net  cooling    (69),    heat  lost  by  long-wave  radiation    (  AQ  ) 
and  total  precipitation    (P)    during  12-hr  period  for  cloud  extending  from 
1000  to  400  mb. 


Initial  6  69  AQ  P 

w  w  9-2 

(deg  C)  (deg  C)  (10     ergs  cm      )         (mm) 


0  1.2  11.2  0.9 

10  0.8  9.8  1.4 

20  0.3  6.1  1.4 


clouds.    The  amounts  increase  slightly  with  initial    0    .     Note  that  a 

w 

cooling  of  only    0.3  C    for    8    =  20  C    produces  more  precipitation  than 

does  a  cooling  of    1.2C      for    0    =0C. 

w 

Precipitation  amounts  given  in  Table  2  are  fairly  modest.    However, 
in  a  quasigeostrophic  numerical  prediction  model,  the  large-scale 
vertical  motion  varies  inversely  as  the  value  assigned  to  the  static 
stability.    Thus  the  cloud  destabilization  may  lead  to  an  increase  in  the 
large-scale  vertical  velocity  also.    This  would  contribute  further  to  the 
amount  of  precipitation. 
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Appendix  -  Derivation  of  Eq  (5) 

Consider  the  equation  for  3-dimensional  motion 

dV  -*       -»      -*      -• 

-=-ttVp-2QxV+g  +  F  (Al) 

where  V   =  u  f  +  v  7  +  wj,    7  =  f  B/dx  +  f  d/Sv  +  k  d/dz   and    F    is  the 
frictional  force  associated  with  motions  arising  from  convective  over- 

-4 

turning.    Multiplying  Eq  (Al)  scalarly  by    V  gives 

£  (  y   +  gz)  =  -  a  V  •  Vp  +  v  •  F  (A2) 

The  first  law  of  thermodynamics  may  be  expressed  as 

r    dT     .  dQ        da  ,     . 

CvTt  *  St    pIi  <*3) 

Adding  (A2)  and  (A3)  and  using  the  equation  of  continuity 

~  =   av  •  v  M) 

yields  the  equation 

2 

d       V  -♦      dO      ->      -* 

—  (y    +gZ  +  CvT)  -fl7'  pV+^+V  -  F  (A5) 

Now  if   F   is  any  function,  use  of  (A4)  gives 

p|f=57  PF  +  V  •  pFV  (A6) 

Integrating  (A5)  with  respect  to  mass  over  the  entire  atmosphere,  and 
making  use  of  (A 6)  and  the  relation 

00  n 

I    p(gZ  +  CT)dz  =  JpCTdz  (A7) 

o  o        p 

one  obtains 
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where    dm  =  pdz  dA  is  an  element  of  mass  . 

Now  if  the  motions  arising  from  convective  overturning  are  dissipated 

as  heat, 

v2 
6  J  -y  dm  =  0  (A9) 

Furthermore , 

j6Qfdm  =  -J[  Jv  •  F  dm]  dt  (A10) 

where    6  Q,  is  the  heat  added  per  unit  mass  due  to  friction  and  the  time 
integration  extends  over  the  interval  to  which    6Q,    refers.    Thus  (A8) 


becomes 


6  J  C  T  dm  =  J6Q  dm  (5) 


where    6  Q    refers  to  heat  added  from  sources  other  than  friction. 
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